Doped t-J Model on a Triangular Lattice: Possible Application to Na^CoC^ • ^H 2 and 
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We report the finding of time-reversal-symmetry-breaking d x 2_ y 2 + id xy superconducting ground 
state in the slave-boson mean-field theory for the t-J model on triangular lattice. For t/J = —5 
(t/J — —9) pairing exists for x < 13% (x < 8%) upon electron doping, and x < 56% (x < 13%) 
upon hole doping. These results are potentially relevant to doped Mott insulators Na^CoC^.yEhO 
and Nai-xTiCb. 
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There is now broad agreement that the physics of 
the high T c cuprate is that of the doped Mott insula- 
tor. However, 17 years after its discovery[l] the layered 
cuprates remain the only materials which exhibit the phe- 
nomenon of high T c superconductivity. There are three 
reasons which make the cuprates unique: (1) The par- 
ent compound is a Mott insulator with S = i and no 
orbital degeneracy; (2) the structure is two dimensional; 
and (3) the exchange energy J is very large (J w 1500 
K). Anderson[2] has stressed the strong quantum fluc- 
tuation of the S — \ system in two dimensions. His 
resonating valence bond (RVB) theory describes a liquid 
of spin singlet which becomes a superconductor when the 
holes are phase coherent. If this is a general property of 
a doped Mott insulator, it will clearly be desirable to 
examine other examples which satisfy these three crite- 
ria. Recently Takada et al.[3] reported the discovery of 
superconductivity in Na^Co02.yH20 with a T c of 5 K 
for x — 0.35. As these authors pointed out, this system 
may be viewed as a Mott insulator with electron doping 
of 35%. The Co atoms are in a triangular lattice and 
the Co 4+ atom is in a low spin (S = |) state. Thus 
this material satisfies the first two criteria listed above. 
The value of J is not known at present, but the new dis- 
covery offers hope that a second system which exhibits 
superconductivity by doping a Mott insulator may be re- 
alized. A summary of what is known about this material 
and a discussion in terms of RVB physics was given by 
Baskaran.[4] In this paper we argue that the t-J model on 
a triangular lattice is a reasonable starting point to model 
these materials. We estimate the value of J, and com- 
pute the slave boson mean field phase diagram for both 
electron and hole doping. We find d x 2_ y 2 + id xy pairing 
over significant range of doping for both t/J = —5 and 
t/J = —9. The appearance of superconductivity shows 
an interesting particlc-holc asymmetry. Wc propose that 
Nai-^TiO^ may be an example of the hole-doped system. 

The transition metal oxide layers in Na^OoC^ and 
Nai-^TiC^ form a common structure where the tran- 
sition metal is surrounded by an octahedral oxygen cage. 
The cages are edge sharing, forming a layered structure. 
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FIG. 1: (a) The NaTi0 2 structure, (b) Projection of Fig.l(a) 
along the (1,1,1) direction showing the oxide layer. This layer 
is similar in Na^OoC^ and NaTi02. Oxygen occupies vertices 
of the solid (dashed) triangles which lie above (below) the 
cobalt or titanium layers, (c) The splitting of the ti g levels 
due to the distortion of the oxygen octahedra. 



The stacking of the oxide layers is different in the two ma- 
terials and we first describe the structure of Nai^TiC^. 
A way to visualize the structure which is convenient for 
understanding the electronic structure is given in Fig. 1(a) 
(adopted from Ref. [5]). It shows a piece of the crystal 
structure in a co-ordinate system which is standard for 
the oxygen octahedron. The hexagonal layered struc- 
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ture is visualized by looking down on this structure along 
the (1,1,1) direction. Thus the axis z' becomes the c- 
axis in the layered hexagonal structure notation. One 
can see stacking of O-Ti-0 or O-Co-0 hexagonal lay- 
ers, highlighted in Fig. 1(b), separated by a Na layer. In 
Nai_ x Ti02, the Na sits directly over the oxygen in the 
layer below Ti (vertices of the dashed triangle), whereas 
in Na a; Co02, the Na site is on top of the Co. The next 
C0O2 layer has Co on top of Na but the stacking order of 
the oxygen layer is reversed (solid triangle now lies below 
Co) thereby doubling the unit cell along c. This differ- 
ence in stacking does not affect the electronic structure 
of the C0O2 layers and will be ignored from now on. 

From Fig. 1(a) it is clear that the d states are split into 
t 2g and e g orbitals by the octahedral environment. A 
shift of the oxygen along the (1,1,1) direction splits the 
t 2g orbitals into [5] 

Alg = (d X y+dy Z +d ZX )/V3 (l) 

and a doublet (labelled E g to distinguish from e g ) 

Eg = [{d ZX - dy Z )/V2, (-2d X y + dy Z + d ZX )/V&) . (2) 

The band calculation of Singh [6] showed that the split- 
ting is fairly large in Nao.sCoC^ with A\ g lying higher 
than Eg. As shown in Fig. 1(c), in Co 4+ the unpaired 
spin occupies the non-degenerate A\ g orbital. We note 
that the t 2g orbitals have lobes which point to the mid- 
point of lines connecting the oxygens in the octahedral 
cage. From Fig. 1(a) and Eq.(l) we see that the A\ g or- 
bitals on nearest neighbor Co have components which 
point directly at each other. Thus unlike the cupratc, 
where the hopping is via the Cu-0 covalent band, in the 
cobalt compounds, the direct overlap between the A\ g 
orbitals form a band. Due to band overlap, it is difficult 
to extract the A\ g bandwidth from Singh's calculations 
[6] but we estimate it to be between 1 to 1.4 eV. From 
this we extract the hopping integral t for the t-J model 




(3) 

to be t = -0.11 to -0.15 cV. In Eq. (3) the pro- 
jection operator P removes double/zero occupancy for 
hole/electron doping respectively. We note that the over- 
lap between A\ g orbitals is positive and the negative 
sign of t is a consequence of the sign convention cho- 
sen in Eq.(3). As emphasized by Baskaran,[4] there is no 
particle-hole symmetry in the triangular lattice and the 
sign of t is important. Singh's band structure [6] shows 
a maximum in the band structure at T, confirming the 
negative sign of t. The Fermi surface of Na .5CoO2 con- 
sists of a hole pocket of area \ around the T point. [6] 
This is consistent with photoemission results. [7] 

It is much more difficult to estimate the exchange con- 
stant J = At 2 /U because the U parameter is highly un- 
certain due to screening. We appeal to another t 2g S = \ 



system where J has been determined experimentally. In 
TiOCl the t 2g orbital is orbitally ordered and forms one- 
dimensional S — \ spin chains. [8] The exchange interac- 
tion was found to be 660 K by fitting the spin suscepti- 
bility to the Bonner-Fisher curve. In this case the t 2g or- 
bitals are pure d yz and point directly at each other. The 
bandwidth from band structure calculations is also about 
1 eV, leading to t w -1/4 cV. Note that t for Na x Co0 2 
is considerably smaller even though the Co-Co distance 
at 2.84 angstrom is shorter than the Ti-Ti distance (3.38 
angstrom). This is because only one component out of 
three in Eq.(l) points directly towards each other. As- 
suming that U is similar, our best guess for J is 12 to 24 
meV and \t/J\ w 6 to 9. Note that our estimate of J for 
Na a; Co02 is about an order of magnitude smaller than 
that for the cuprates. 

It is interesting to consider another S = | system 
which corresponds to single occupation of the t 2g orbitals. 
NaTi02 has similar layer structure as Naa;Co02 but the 
Na layer is nominally fully occupied. The Ti 3+ has d 1 
configuration giving rise to unpaired S = \. Contrary 
to Nao.5Co02, the Ai g , E g splitting is much smaller than 
the bandwidth, according to LDA calculations. [5] How- 
ever, an LDA+U calculation shows that the A\ g orbital is 
occupied preferentially, so that the single electron again 
occupies the nondegenerate A\ g band. [5] Reduction of 
the Na occupation (Nai_2,Ti02) corresponds to doping 
by x holes. The undoped system is of great interest be- 
cause it is one of the few known examples of the S = 5 
triangular antiferromagnct. However, the control of Na 
stoichiometry presents a serious materials challenge and 
relatively few studies have been carried out to date. [9] 

As we mentioned earlier the occupation constraint of 
Eq. (3) is different for the electron and hole doping. In 
the case of electron doping we shall perform a particle- 
hole transformation Cj l(7 — > cj _ CT so that the constraint 
always means no double occupancy and the sign of t is 
reversed. Consequently t < for hole-doped Nai_. E Ti02 
and t > for electron-doped Naa;Co02. 

The starting point of our calculation is the following 
U(l) slave-Boson mean-field Hamiltonian[10] for the t-J 
model at hand, H MF = H nm + H m with 

(ij)crij' i<y 

H m =J^{S i -m j + S j -m i ). (4) 

(ij) 

Here xi — 1 — 'Yl ia {ni a ) and fj,i is the correspond- 
ing local Lagrange multiplier. The doping level is 
given by x = J2i x i/N with N the lattice sites. The 
mean field order parameters are Xij — ^2 a (ftofj<r)> 

= Xw W(/i<7 •/>'>! and m i = ( S i) with S i = 
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FIG. 2: Dependence of the order parameters on hole doping 
(a;): t = 5 J (a), 9 J (b), -5 J (c), and -9 J (d). Note the scale 
change in 2(b) and 2(c). 



J2 SS > <Jss ' flsfis' /2- Note that a general non-collinear 
magnetic order is allowed. 

Our search for the mean-field solution proceeds in two 
stages. First we assume the density to be uniform and 
perform an unrestricted search of the parameters Xij > Ay 
and m.i (over 6x6 elementary unit cells in a large lattice). 
For x < x c we find Xij to be a real constant as long as x 
is nonzero and Ay = while the spins form the V3 x \/3 
non-collinear structure identified by Husc and Elser[ll] 
for x = 0. If we denote the lattice by = na + Ih, 
a = x and b = — ^x + ^y, then a representative so- 
lution has mf = and mf + im? — |m|e lC *' ri , where 
Q = (27r/3,27r/V3) and Q r, = ^(n + I). |m| = § 
at zero doping, and decreases gradually with increas- 
ing doping. It is roughly 0.35 ~ 0.38 (depending on 
our choices of t/J) at x — x c , where the long-range 
ordered magnetic state undergoes a first order phase 
transition into a uniform d x 2_ y 2 ± id xy superconduct- 
ing state. The value of x c is 3.8%, 3.7%, 2.2%, 2.1% for 
t/J = 5,— 5,9,— 9 respectively. Even if the pairing chan- 
nel is removed, the magnetic order extends only up to 
Xt = 4.5%, 4.2%, 2.4%, 2.4% for the same set of t/J. 
This is significantly different from the mean-field solu- 
tion of the square lattice with t/ J = 3 as appropriate for 
the cuprates. In that case the system is inhomogeneous 
with alternating stripes of superconducting and magnetic 
regions at roughly 2% < x < 12%, and the magnetic or- 
der becomes commensurate with coexisting uniform pair- 
ing order at higher doping levels up to x c <~ 20%. [12] 
For the hexagonal lattice, the stable region of the anti- 
ferromagnetic solution is greatly reduced. This is rea- 
sonable in view of the frustration of the lattice, which 
disfavors magnetic order and was the original motivation 
of the RVB as a competing state [2]. 



In the second stage of the calculation we concentrate 
on doping greater than x c . In view of the first-order na- 
ture of the transition to magnetic order, we shall present 
results without magnetic order even for x < x c , with the 
understanding that the solution is only locally but not 
globally stable. We suppress magnetic order and per- 
form an unrestricted search for solutions of H nm in Eq. 
(4) on a 400 x 400 lattice. In the elementary unit cell 
(varied from 1 x 1 up to 4 x 8) we put random initial 
values for Ay, x i an d Mi, and evolve them so that 
the mean-field free energy is minimized at the desired 
doping level. Note that this time we allow nonuniform 
charge density. However, with our choice of parameters, 
we found spatially non- uniform solutions only for x < x u , 
where x u is always less than x c and these will be ignored 
from now on. Due to the uncertainty of t/J we have 
studied two cases \t/J\ = 5 and \t/J\ = 9. It turns out 
that in both cases there is a significant doping range in 
which the ground state is superconducting. The pairing 
symmetry is always d x 2_ y 2 + id xyi hence breaks the time 
reversal symmetry. It is interesting to note that if one 
ignores the magnetic order parameter at zero doping the 
same mean-field theory predicts a degenerate family of 
ground states. Among them the d x 2_ y 2 + id xy paired 
state and the 7r/2 flux state are two examples. [13, 14] To 
visualize the pairing pattern, consider an arbitrary site. 
The associated with the six bonds stemming from it 



has the form A, 



,j20, 



for d x 2_ y 2 + id xy pairing and 



Ay 



Ae _t29i3 for d x 2_ y 2 — id xy pairing. To generate the 
pairing field over the entire lattice, simply translate the 
above pattern to other sites. In the absence of a mag- 
netic field the d x 2_ y 2 +id xy and d x 2_ y 2 — id xy pairing are 
degenerate. 

The doping dependence of the zero temperature pair- 
ing order parameter A = |Ay| is plotted in Figs. 2 (solid 
lines) for t = 5J (a), 9J (b), -5J (c), and -9J (d), 
where we observe that pairing exists for x < 13% (a), 8% 
(b), 56% (c) and 13% (d). The peak in A at x = 0.5 in 
Fig. 2(c) is explained by the Van Hove peak in the free 
electron density of states. Also shown in Figs. 2 is the 
hopping order parameter \ — I Xij I (dashed lines). It is 
weakly dependent on x for x < 30%. Interestingly Xij 
has the same sign of t from our calculation. This is a 
reasonable result as it increases the band width of the 
fermions so that the kinetic energy is lowered. Finally, 
a much more important feature in Figs. 2 is the consid- 
erable asymmetry between electron doping (t > 0) and 
hole doping (t < 0). This is due to the particle- hole 
asymmetry in the free-electron dispersion on the trian- 
gular lattice. Indeed, the stronger pairing for hole doping 
is due to the increases in the Fermi level density of state 
as the averaged occupation decreases from half-filling. 

We have also computed the onset temperature of the 
RVB fcrmion pairing and Bosc-Einstcin condensation 
(BEC) of holons.[15] The mean-field superconducting 
transition is the smaller of the RVB and BEC curves 
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FIG. 3: Dependence of the RVB and BEC critical tempera- 
tures (in units of J) on hole doping (x): t = 5 J (a), 9 J (b), 
— 5J (c), and — 9J (d). Note the scale change in 2(b) and 
2(c). 



(Figs. 3). As usual, the Bose condensation temperature 
is an overestimate of the phase coherence temperature of 
the slave bosons. More generally because of the prox- 
imity to the Mott insulator limit, we expect that the 
superconducting transition temperature will be deter- 
mined by the supcrfluid density at low doping, and by 
the onset of pairing at higher doping, forming a phase 
diagram similar to that of the cuprates. It is worth not- 
ing that in the case of cuprates, the thermal excitation 
of nodal quasi-particles significantly reduces the super- 
fluid density and therefore the transition temperature at 
low doping. [16, 17] The existence of a full gap in the 
d x 2_ y 2 +id xy state suppresses this possibility and every- 
thing else being equal, we can expect a higher T c in this 
case. 

We conclude that within the slave boson mean field 
theory, the systems exhibits non-colincar antifcrromag- 
nctism for x < x c . Above this lower critical concentration 
and below certain upper critical concentration the system 
exhibit time reversal symmetry breaking d x i_ y i + id xy 
pairing state. The orbital moment produces a magnetic 
field which can be detected by /iSR. The field has been 
estimated to be 15 Gauss for the anyon model [18] but 
a calculation based on the slave boson mean field the- 
ory yields a smaller estimate of 1 Gauss near a vacancy. 
[19] If the orbital moment in neighboring layers are par- 
allel, this state has been predicted to exhibit fascinating 
new effects such as quantized spin Hall conductance and 
anomalous Hall thermal conductivity. [20] The orbital mo- 
ment corresponds to roughly 4j^-B, and field cooling in a 
modest magnetic field may line up the orbital moments. 
It will clearly be desirable to achieve such a state in the 
laboratory. For \t/J\ — 5 (\t/J\ — 9), the d x i_ y i + id xy 
superconductor is expected to exist over a low doping 



range x < 13% (x < 8%) for electron doping, and a wider 
range x < 56% (x < 13%) for hole doping. Thus some 
ingredient in addition to the t-J model may be needed to 
explain the experimental observation of superconductiv- 
ity at x — 0.35 in electron-doped Na x Co02-j/H 2 0. Explo- 
ration of the Nai-zTiO^ system may be more promising. 
We also note that Nao.7Co02 is a metal which exhibits 
unusual behavior such as linear T resistivity and large 
magnetic-field-dependent thermal power. [21] This is also 
inconsistent with the mean field prediction of Fermi liq- 
uid in the overdoped region, suggesting that some ad- 
ditional physics may be at work. In any event, the new 
observation opens up the possibility of changing the dop- 
ing concentration in a controlled way by intercalation and 
much new physics surely remains to be discovered. 

After the completion of this work we have seen a paper 
by Kumar and Shastry [22] reaching similar conclusions. 
However, we do not agree with their identification of the 
sign of t. 
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